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At low temperatures, the transport through a superconducting-normal tunnel 
interface is due to tunneling of electrons in pairs. The probability for this process 
is shown to depend on the layout of the electrodes near the tunnel junction, rather 
than on properties of the tunnel barrier. This dependence is due to interference of 
the electron waves on a space scale determined by the coherence length, either in the 
normal or the superconducting metal. The approach developed allows us to evaluate 
the subgap current for different layouts of interest. 
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It is well-known that the charge transport through a tunnel NS interface between a 
normal metal and a superconductor is strongly suppressed at voltages lower than A/e, A 
being the superconducting energy gap [[!]. Indeed, energy conservation forbids the transfer 
of a normal electron with an energy below the gap to the superconductor, since it would 
have been converted into a quasiparticle with an energy larger than A. 

Experimentally, some residual conductivity has been observed at subgap voltages even 
at very low temperatures. There is a tendency to ascribe this either to imperfections in the 
tunnel barrier or to normal inclusions in the superconductor. Another mechanism of the 
subgap conductivity is the so-called two-electron tunneling @. The point is that two normal 
electrons can be converted into a Cooper pair, thus this transfer may cost no energy. The 
current will be proportional to the fourth power of tunnel matrix elements; therefore it is 
much smaller than the one-electron current. 

The problem was previously treated under the simple assumption that the electron wave- 
functions in both metals are just plane waves. In this case one can consider a barrier of 
arbitrary transparency in order to describe the crossover from a tunnel to a perfectly con- 
ducting interface [[J. But some important physics may be missed under this assumption. 
Let us compare the two realizations of the NS interface depicted in Fig. 1. In Fig. la the 
electron transmitted through the interface does not experience any scattering in the metals 
and never gets back to the junction. The plane wave picture seems to be relevant for such a 
geometry. An alternative situation is shown in Fig. lb. This case is usually realized when 
the tunnel junction is formed by imposing two thin metal films. The transmitted electron 
gets back to the junction many times before leaving the junction region. Thus the tunneling 
occurs between electron states of very complex structure which emerges from interference 
between scattered waves. 

This interference has no effect on one-electron transport, since the average one-electron 
density of states does not depend on the stucture of the wave function. However, it matters 
for two-electron tunneling, since two electrons penetrating the barrier will interfere. Such an 
interference occurs at a space scale corresponding to the energy difference between the two 
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electron states. It makes the probability of two-electron tunneling dependent on the system 
layout at the corresponding mesoscopic space scale. Our aim is to evaluate this interference 
effect for an arbitrary given layout. The rapid progress of nanotechnology makes it possible 
to fabricate numerous relevant stuctures, so it is worthwile to be able to give guidelines 
to a designer. As we will see below, the subgap conductivity is strongly enhanced if the 
interference effect is essential. 

We first review shortly the two-electron tunneling through a superconducting-normal 
interface as it has been discussed by Wilkins and more recently by Hekking et al. [|J. 
The total Hamiltonian can be written as H = H^+Hs+Ht- The subscripts N and S refer to 
the normal and the superconducting electrode respectively; the transfer of electrons through 
the tunnel interface is described by the tunnel Hamiltonian H?. The latter is expressed in 
terms of quasiparticle operators 7, 7^ for the superconductor, and electron operators a, a* 
for the normal metal: 

Ht = {*k P ak, ff (Wp, CT 7p,<7 + v Pja jl p _ a ) 

k,p,cr 

+*kpOW7 P>(7 + f P ,«r7-p,- ff )ak,4- (1) 

Here, tk p are the tunnel matrix elements which we take to be spin-independent, and w P)CT , f P)(T 
are the BCS coherence factors [p]]; the sum is taken over momenta k, p and spin a =f, \. 

Using second order perturbation theory in Ht one can calculate the amplitude for the 
transfer of two electrons from the normal to the superconducting electrode: 

^kTk'i = E Vk'-p^P I J^— + T^— ) • (2) 

Here the spin dependence of the coherence factors was dropped after using the relation 
f p l = — w_ p ,j.. We define electron energies and ( p for the normal and the superconducting 
electrode respectively, and quasiparticle energies e p = JA? + Q. The denominators in (0) 
reflect the fact that a virtual state is formed when the first electron enters the superconductor 
as a quasiparticle. The second electron couples to this quasiparticle, thus forming a Cooper 
pair. The corresponding rate r(V) as a function of the voltage V applied across the junction 
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can be found by using Fermi's Golden Rule 

4<7T 

T ( v ) = IT E \ A Wl\ 2 f(Zk)f(Zk>m k + & + 2eV). (3) 
n k,k' 

It contains the Fermi functions / for electrons with energies in the normal metal. We 

recall that the normal conductance of the junction is determined by the rate 7(V) for usual 
electron tunneling, which is proportional to |£k P | 2 : j{V) — |£k P | 2 ,/(£fc)(l — /r(Cp))^(£ft — Cp + 

e y). 

The calculation of | ^4ktk'| | 2 i n Eq. @ involves summations over momentum of a product 
of four tunnel matrixelements. It therefore requires an assumption about the dependence 
of the tkp on the wave vectors k and p. This dependence is strongly related to the nature 
of the electron motion in the electrodes, as we discussed above. Following || we assume 
first that plane electron waves propagating in the electrodes are transmitted specularly by 
a rectangular tunnel barrier with, say, a length Lj, and a height U (See Fig. la). The area 
of the junction will be denoted by S. Specular scattering implies that the components of 
momentum ky and pii parallel to the barrier plane are conserved. If S is of the order of A^ 
(with Xp the Fermi wavelength of the electrons) the values of ky and py are quantized, leading 
to discrete transport channels ||. The corresponding quantum numbers are equal: = 
n p ; the effective number of channels N e ff contributing to the transport will be calculated 
below. The magnitude of tk P decreases exponentially with decreasing squared component 
k^ perpendicular to the barrier. This results in the assumption t kp oc <5 nki „ p exp — k|A 2 , 
with A = hLb/y/8mU, where m is the electron mass. The calculation of the rate @ is 
easily performed using this model for tkp, by averaging products of these matrix elements 
over directions of momentum. As a result we find that V oc Similarly we obtain 

7 oc 1/N e ff. Comparing T with j 2 we find the effective number of transport channels 
penetrating the barrier, N e ff. 

£ oc (|(tkptk' P > p | 2 W (1^1% = 4ttA 2 /S = l/N eff . 

This result is obtained by assuming ballistic motion of the electrons in the electrodes. 
This assumption is correct only if the scattering of the electron is negligible. Scattering 
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may occur at the boundaries of the electrodes or at impurities inside the electrodes. Both 
processes can be characterized by a space scale l e , which corresponds to the distance the 
electron traverses before undergoing the first scattering event. Interference occurs on a 
space scale £ cor . The ballistic picture is valid if the typical size of the junction or £ cor is 
smaller than l e . When these lengths are of the same order, we expect a cross-over to different 
behaviour. In the opposite limit the electron moves diffusively in the junction region. Due to 
interference between incoming and backscattered electron waves N e ff will decrease, thereby 
increasing the rate T, and hence the conductance due to two-electron tunneling. 

Now we will present a method to describe two-electron tunneling in the diffusive transport 
regime, employing the quasiclassical approximation. This enables us to evaluate the tunnel 
matrix elements and express the rate T in terms of quasiclassical diffusion propagators. 
The method is similar to the one presented in Ref. j7|. We start by rewriting the matrix 
elements tkp = / drdr'ipl(r)ij)p(r')t(r,r'), where ipp{ r ) forms a complete set of one-electron 
wave functions in the electrodes, and t(r, r') describes the tunneling from a point r' in 
the superconductor to a point r in the normal metal (primed space arguments refer to 
the superconductor). We also define a propagator from r 2 to r\ by K^(ri,r 2 ) = Y,k$(£ ~ 
6c)V , fc( r i)V ; fc( r 2)- With these definitions it is possible to rewrite Eq. @ as 
4tt r 

iW = T J dtdZ'd(d('F((; £ one'; £, Os(C, C; & O/(0/(O*(£ + + ^) (4) 

with F(C; £, e) = u(C)v(C) {(^ - e)- 1 + ({' - s)- 1 } where e = V^+C, and 

S(C,C';£,0 = / d 3 n...d 3 n J dM-^M**(n,ri)t*(r 2 ,^)t(r 3 ,^)t(r4,ri)x 

^(r 1 ,r 3 )^(r 2 ,r4)K c (r^ri)Xc'(^3 ) ^) (5) 

The physical meaning of Eq. (^) can be understood easily by depicting the integrand of 
Eq. d^D diagrammatically, as has been done in Fig. 2. We see two electrons that propagate 
in the normal electrode with energy £ and The first electron reaches the barrier at ri 
and tunnels to r[, the second electron tunnels from r 2 to r 2 ; both change their energy to 
(. In the superconductor they form a Cooper pair. Since tunneling occurs only between 
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neighboring positions, we have in addition T{ ~ r-. The diagram expresses a probability, and 
therefore is completed by adding the time-reversed process. 

To analyze expression ([5]), it is important to consider the scale of separation of the 
coordinates r\ ~ r[, r 4 ~ r 4 lying on the interface. In the ballistic transport regime, these 
coordinates are separated only by a few Fermi wavelengths. In this case the contribution 
depends on properties of the tunnel barrier only. Below we will concentrate on contributions 
to (HI) which arise when the region of integration is defined by coordinates which are pairwise 
close, but with the pairs separated by a distance much larger than the Fermi wavelength. 
These contributions contain averaged products of two propagators K, which are known to 
decay on a mesoscopic scale in the diffusive transport regime . These products correspond 
to the semiclassical motion of electrons from one point on the interface back to another point 
on this interface. They describe the interference between scattered waves. In the diffusive 
regime these contributions dominate; that is why we concentrate on them. 

There are three types of contributions, which are depicted in Fig 3. Fig. 3a corresponds 
to the case r\ ~ r 2 and r 3 m r 4 . The interference contribution originates from the normal 
electrode. Fig. 3b describes the opposite situation with interference occurring in the super- 
conducting electrode. Here, T\ « r 3 and r 2 ~ r 4 . Finally, in Fig. 3c, interference occurs both 
in the normal and in the superconducting electrode, when t\ ~ r 4 and r 2 ~ r 3 . However, 
we estimated this contribution to be less important, compared to the previous diagrams. 
Therefore, the total effect can be represented as the sum of the interference contributions 
from the superconducting and the normal metal. 

As an example we will discuss the contribution of Fig. 3b to the rate (|]) in some detail. 
The averaged product of the propagators in the superconductor determines the semiclassical 
conditional probability P(r[,r' 2 ]n,n']t) that an electron with position r' 2 and momentum 
direction n' at time t — has position r[ and momentum direction n at time t. Since the 
tunnel amplitude t(r,r') is nonzero only when r and r' are close to the junction interface, 
we can restrict spatial integrations to planar integrations over the junction surface. It is 
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possible to show that 

~ s = ~ 5 (C - CO = tt^j- I d 2 r[d 2 r> 2 f d 2 nd 2 n' g{n,r' x )g{n' ,r' 2 )x 

J dte^-^P^r'^n'-t) (6) 

where vs is the density of states for the superconductor for two spin directions and / d 2 n 
denotes integration over directions of momentum. The function g(n,r) defines the normal 
conductance of the junction: Gt = J d 2 r J d 2 ng(n,r). An expression similar to (^) can be 
obtained for Fig. 3a, by replacing subscript S — > N, energies C — *■ £, and primed space 
arguments by unprimed ones. 

Let us start our consideration of concrete layouts with the simplest geometry of an infinite 
uniform junction between a normal and a superconducting film (Fig. 4a). We assume that 
the film thickness is much less than the coherence length in the superconductor. Then we 
can exploit the picture of two-dimensional electron diffusion. The probability function we 
need is given by 

d being the thickness of either the superconducting or the normal metal film. Taking the 
Fourier transform of this function with respect to time and integrating (Q) over coordinates 
we obtain 

2s(C - CO = (^e'Susds)-^ - CO; Zn(Z - CO = {^ 2 e'Sv N d N )-H^ - CO; (8) 

The current is given by a sum of two terms (eV ~ A ^> T): 

I(V) = I N + I s ; 
2G 2 T h 

In 



e 3 Sis N d 



N 

Glh eV 

*s = / ( 9 ) 

rSusdsnAiJi-eV/A 

It is plotted in Fig. 5. The part emerging from the interference in the normal metal does 
not depend on voltage. So the current sharply jumps at zero voltage, provided T = 0. The 
jump is smoothened at voltages of the order of the temperature: 



I(V, T) = I N tanh(eV/2T) (10) 

The other contribution diverges near the threshold voltage indicating the necessity to make 
use of higher order terms in tunneling amplitudes to describe the crossover between two- 
electron and one-electron tunneling. 

It is worthwile to compare the magnitude of the result with the one we derived assuming 
ballistic motion. The order of the ratio at voltages of the order of A is Imt/haii — idean/d 
i Mclean being the coherence length in the pure superconductor. Therefore the interference 
term dominates under the assumptions we made. 

The coherence along a normal or superconducting film is characterized by the coherence 



length = ^D/eV, JD/A, respectively. The relations @ are valid, provided the junc- 
tion size is much larger than these lengths. In the opposite limit of small junctions, the 
subgap conductivity will be determined by the junction surroundings, rather than by the 
junction itself. Let us illustrate this by considering the geometry in Fig. 4b, where a normal 
electrode is connected to a superconducting sheet by the tunnel junction. In this case we 
find 

S*(C - CO = In TT—WT I S*(* - O = In 7^7^ (n) 

e z 87T 4 ((, — Q')t e A 87r 4 (£ — £')r 

The time r is of the order of S/D, the time spent in the junction area, and provides a 
lower cut-off for the fourier integral. The sheet resistance of the normal (superconducting) 
film is given by = (e 2 zAD7v(s , )^7V(s , ))~ 1 - Indeed does the result not only depend on the 

properties of the tunnel barrier itself (through the dependence on Gt) , but also on properties 
of its surroundings through the dependence on R n . Moreover, the dependence on the precise 
geometry of the layout enters through numerical prefactors. If, e.g., the tunneling would 
occur towards an infinite superconducting sheet instead of a superconducting half plane, the 
semiclassical probability P would be twice smaller, thus decreasing S, and hence the rate 
T, by a factor of 2. The current is again given by a sum of two terms In and Is' 

In = — Rn,NG?p\n——- ; I s = — Rn tS G^\n—— (12) 

7T eVT 7T Ar 
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Note that, in contrast to Eq. (|9|), the subgap conducticity depends only weakly on the 
junction area through the cut-off time r. 

We finally consider the geometry depicted in Fig. 4c. It consists of a small island (length 
Ls, thickness d,s), coupled to two macroscopic leads by tunnel barriers. The grain is linked 
capacitively to the leads. Electron transport through such a system, characterized by a 
small electric capacitance C, has been studied both experimentally and theoretically in great 
detail during the past years f|. The key point is that variations of the charge of the island 
in the course of electron tunneling increase the electrostatic energy, typically by an amount 
Ec = e 2 /2C. This is why electron tunneling through a small grain is suppressed (Coulomb 
blockade). The case of a superconducting island connected to two normal electrodes (NSN 
geometry), was studied recently in Refs. [p[,|9UlO| . Our method to include interference effects 
can also be applied to this case. The charging energy Ec will enter our results explicitly via 
the virtual state denominators in (|) j|, resulting in a dependence of the function F on Ec- 
We will restrict ourselves to the situation in which Ec is smaller than the superconducting 
gap: Ec < A. In order to calculate the contribution due to interference on the island we 
assume that the time h/(A — Ec) spent by the virtual electron on the island is much longer 
than than the classical diffusion time L 2 S /D. If A < Ec, the size of the island is smaller 
than ^ or . In this case, the electron motion covers the whole island and the probability P is 
constant: P = l/(L|dg). As a result we find: 

where ws = l/{vsL 2 s ds) denotes the level spacing of the island, which shows once more that 
the rate (|4]) is not only determined by properties of the tunnel barrier. The corresponding 
current Is reads 

/s = V « G |^ x 
ixe 2 E c 



7T 2A / El \ IA-Ec , AE C 
1 1 arctan W — — — h 



2 



2 JA^ - El I ^ 2 -El\ ^A + E c A*-E c 



(14) 
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When interference in the normal electrode is taken into account, we find Hjv(£ — £') = 
(h/e 2 )(R%G 2 T /4n 4 ) ln[n/(£ - f'H lik e in Eq. (Jig), however larger by a factor of 2, due the 
difference in geometry. The resulting current reads: 

I N = %R N n Gl In 4- l 1 ^ arctan (15) 

In conclusion, we evaluated the low-voltage supgap conductivity of NS boundary inter- 
face. In many interesting cases it was shown to be determined by the conditions of electron 
motion in the electrodes rather than the properties of the tunnel barrier. Therefore it de- 
pends on the system layout on the mesoscopic scale. We presented an approach which gives 
exact results for any layout given. 

The authors are grateful to M. Devoret, D. Esteve, J. Mooij, A.Schmid, and G. Schon 
for a set of valuable discussions. 



10 



REFERENCES 

[1] M. Tinkham, Introduction to Superconductivity (McGraw-Hill, New York, 1975). 

[2] J.W. Wilkins, in Tunneling Phenomena in Solids, edited by E. Burstein and S. Lundqvist 
(Plenum Press New York, 1969), p. 333. 

[3] G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. Rev. B 25, 4515 (1982). 

[4] F.W.J. Hekking, L.I. Glazman, K.A. Matveev, and R.I. Shekter, submitted to Phys. 
Rev. Lett. 

[5] Y. Imry in Directions in Condensed Matter Physics, edited by G. Grinstein and G. 
Mazenko (World Scientific, Singapore, 1986), p. 100 

[6] B.L. Al'tshuler and A.G. Aronov in Electron-electron interactions in disordered systems, 
edited by A.L. Efros and M. Pollak (North-Holland, Amsterdam, 1985), p. 21 

[7] D.V. Averin and Yu.V. Nazarov in Single Charge Tunneling, edited by H. Grabert and 
M.H. Devoret (Plenum Press, New York, 1992), Chapter 6.3 

[8] D. V. Averin and K. K. Likharev, in Mesoscopic Phenomena in Solids, edited by B.L. 
Altshuler et al. (Elsevier, Amsterdam, 1991), p. 173; Single Charge Tunneling, edited 
by H. Grabert and M. H. Devoret (Plenum, New York, 1992). 

[9] M. T. Tuominen, J. M. Hergenrother, T. S. Tighe, and M. Tinkham, Phys. Rev. Lett 
69, 1997 (1992). 

[10] T. M. Eiles, J. M. Martinis, and M. H. Devoret, submitted to Phys. Rev. Lett. 



11 



FIGURES 

FIG. 1. Two typical realizations of a tunnel junction between a normal (N) and a supercon- 
ducting (S) electrode. In (a) the electron moves ballistically, in (b) diffusively in the junction 
region. 

FIG. 2. Diagram corresponding to Eq. (|5|). Electrons propagate (solid lines) with energy 
in the normal (N) and energy £' in the superconducting (S) electrode. They tunnel through the 
barrier (shaded region) at postitions 1, ...,4, marked by crosses. 

FIG. 3. Contributions to the subgap conductivity due to interference in (a) the normal elec- 
trode, (b) the superconducting electrode, and (c) both electrodes. 

FIG. 4. Various layouts discussed in the text: (a) infinite uniform junction, (b) normal electrode 
connected to a superconducting halfplane, and (c) a superconducting island connected to two 
normal electrodes. 

FIG. 5. I — U-characteristcs for an infinite junction. The curves (from bottom to top) represent 
I s , In, and I(V) = I s + I N 
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